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Abstract. We study the relaxation of a test particle immersed in a bath of field particles interacting via 
weak long-range forces. To order 1/N in the N — > +00 limit, the velocity distribution of the test particle 
satisfies a Fokker-Planck equation whose form is related to the Landau and Lenard-Balescu equations in 
plasma physics. We provide explicit expressions for the diffusion coefficient and friction force in the case 
where the velocity distribution of the field particles is isotropic. We consider (i) various dimensions of space 
d — 3,2 and 1 (ii) a discret spectrum of masses among the particles (iii) different distributions of the bath 
including the Maxwell distribution of statistical equilibrium (thermal bath) and the step function (water 
bag). Specific applications are given for self-gravitating systems in three dimensions, Coulombian systems 
in two dimensions and for the HMF model in one dimension. 

PACS. 05.20.-y Classical statistical mechanics - 05.45.-a Nonlinear dynamics 



1 Introduction 

Kinetic theories of many-particles systems are important 
to understand the dynamical evolution of the system and 
to determine transport properties. The first kinetic equa- 
tion for a Hamiltonian A^-body system was derived by 
Boltzmann in his theory of gases P^. In that case, the 
particles do not interact except during strong collisions. 
Kinetic theories were later extended to the case of parti- 
cles in interaction by Landau [2] in the case of plasmas 
and by Chandrasekhar |3] in the case of stellar systems. 

' In this Introduction, we present a short historical review 
of kinetic equations with gravitational or coulombian in- 
teraction, stressing the literature before the sixties when 
these topics were developed. Additional references can be 
found in the books of Balescu |3j and Saslaw jS], or in 

\ the review of Kandrup |S] . The rest of the paper considers 
extensions of these kinetic theories to new situations. 

Landau |2] derived his kinetic equation by starting 
from the Boltzmann equation and considering a weak de- 
flection limit. Indeed, for a Coulombian potential of inter- 
action slowly decreasing with the distance as r -1 , weak 
collisions are the most frequent ones. Each encounter in- 
duces a small change in the velocity of a particle but the 
cumulated effect of these encounters leads to a macro- 
scopic process of diffusion in velocity space. This treat- 
ment, which assumes that the particles follow linear tra- 
jectories with constant velocity in a first approximation, 
yields a logarithmic divergence of the diffusion coefficient 
for both small and large impact parameters but the equa- 
tion can still be used successfully if appropriate cut-offs 



are introduced. A natural lower cut-off, which is called 
the Landau length, corresponds to the impact parame- 
ter leading to a deflection at 90°. On the other hand, in a 
neutral plasma, the potential is screened on a distance cor- 
responding to the Debye length. Phenomenologically, the 
Debye length provides an upper cut-off. Later on, Lenard 
7 and Balescu |S] developed a more precise (but also 
more formal) kinetic theory that could take into account 
collective effects. This gives rise to the inclusion of the 
dielectric function e(k,k ■ v)| 2 in the denominator of the 
kinetic equation. Physically, this means that the particles 
are "dressed" by a polarization cloud. The original Lan- 
dau equation, which ignores collective effects, is recovered 
from the Lenard-Balescu equation when |e(k, k ■ v)| 2 = 1. 
However, with this additional term, it is found that the 
logarithmic divergence at large scales is now removed and 
that the Debye length is indeed the natural upper length- 
scale to consider. 

In stellar dynamics, Chandrasekhar |3] developed a ki- 
netic theory in order to determine the timescale of col- 
lisional relaxation. He computed in particular the coeffi- 
cients of diffusion and friction (second and first moments 
of the velocity increments) by considering the mean ef- 
fect of a succession of two-body encounters. Since this ap- 
proach can take into account large deflections, there is 
no divergence at small impact parameters and the grav- 
itational analogue of the Landau length appears natu- 
rally. However, this approach leads to a logarithmic di- 
vergence at large scales that is more difficult to remove 
than in plasma physics because of the absence of De- 
bye shielding for the gravitational force. In a series of 
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papers, Chandrasekhar & von Neumann [""] developed a 
completely stochastic formalism of gravitational fluctua- 
tions and showed that the fluctuations of the gravitational 
force are given by the Holtzmark distribution (a particular 
Levy law) in which the nearest neighbor plays a prevalent 
role. From these results, they argued that the logarithmic 
divergence has to be cut-off at the interparticle distance. 
However, since the interparticle distance is smaller than 
the Debye length, the same arguments should also apply in 
plasma physics, which is not the case. Therefore, the con- 
clusions of Chandrasekhar & von Neumann [""] are usually 
taken with circumspection. It is usually argued (e.g., Co- 
hen et al. [ID], Saslaw [Hj, de Vega & Sanchez [I"Q) that 
the logarithmic divergence should be cut-off at the physi- 
cal size R of the cluster for finite systems or at the Jeans 
scale for infinite systems, since the Jeans length is the pre- 
sumable analogue of the Debye length in the present con- 
text (see Kandrup [""])■ Chandrasekhar [12] also developed 
a Brownian theory of stellar dynamics and showed that, 
on a qualitative point of view, the results of kinetic the- 
ory could be understood very simply in that framework. 
In particular, he showed that the coefficients of diffusion 
and friction are related to each other by an Einstein re- 
lation ["d]. Later on, Rosenbluth et al. [*""""| proposed a 
simplified derivation of the coefficients of diffusion and 
friction for plasmas and stellar systems and, substituting 
these expressions in the general form of the Fokker-Planck 
equation, they obtained a kinetic equation which is at the 
basis of the dynamics of stellar systems. They also pro- 
vided simplified expressions of this equation in the case of 
axial symmetry. 

It is interesting to note that the previous authors did 
not point out the link with the Landau equation and 
that modern textbooks of astrophysics [IH] usually derive 
the kinetic equation of stellar dynamics from the Fokker- 
Planck equation by using the approach of Chandrasekhar 
PI and Rosenbluth et al. [Il|. We will see that we can 
cquivalcntly obtain the kinetic equation of stellar dynam- 
ics from the Landau equation. This alternative derivation 
can re-inforce the connection between stellar systems and 
plasmas. We note, however, an important difference be- 
tween stellar dynamics and plasma physics. Neutral plas- 
mas are usually spatially homogeneous due to the Debye 
shielding. By contrast, stellar systems are inhomogeneous. 
Therefore, the above-mentioned kinetic theories developed 
in astrophysics rely on a local approximation. The collision 
term is calculated as if the system were homogeneous or 
as if the collisions could be treated as local. Then, the 
effect of inhomogeneity is taken into account in the ki- 
netic equation by introducing an advective term (Vlasov 
term) in the left hand side which descibes the evolution 
of the system due to mean-field effects. The local approx- 
imation is supported by the stochastic approach of Chan- 
drasekhar & Von Neumann ["J] showing the preponderence 
of the nearest neighbor. However, this remains a simplify- 
ing assumption which is not easily controllable. It is likely 
that the logarithmic divergence at large scales comes from 
this approximation. More recently, Kandrup [TS] derived a 
generalized Landau equation by using projection operator 



technics. This formal approach is interesting because it 
can take into account effects of spatial inhomogeneity and 
memory which are neglected in the previous approaches 
. It clearly shows which approximations are needed in 
order to recover the Landau equation. However, the gen- 
eralized Landau equation remains extremely complicated 
for practical purposes. 

Until now, the kinetic theories of systems with long- 
range interactions have been essentially developed for 3D 
systems with Coulombian or Newtonian potential. The 
main object of this paper is to extend these theories to 
other dimensions of space d — 2 and d = I and for a wide 
class of potentials of interaction. This generalization has 
been initiated in Chavanis |21I22| . It is shown that for sys- 
tems with weak long-range interactions, the Landau and 
Lenard-Balescu equations describe the collisional dynam- 
ics of the system to order I/TV in a proper thermodynamic 
limit N — > +oo. Therefore, the collisional relaxation time 
scales in general as tu ~ Ntjj, where to is the dynamical 
time. For N — > +oo or t <C tu, the collision term is negli- 
gible and we obtain the (mean-field) Vlasov equation. In 
Rcf. 22 , we have introduced a Fokker-Planck equation 
which describes the evolution of a test particle in a bath 
of field particles. We have obtained analytical expressions 
of the diffusion coefficent and friction force in the case of a 
Maxwellian distribution of field particles (thermal bath) . 
We shall obtain here generalizations of these results. 

In the case of a Newtonian or Coulombian potential 
of interaction in d = 2 the diffusion coefficient diverges 
linearly at large scales so that an upper cut-off (related 
to the Debye length in the case of a plasma) must be in- 
troduced. We have suggested in [22] that this divergence 
would be removed if we use the Lenard-Balescu kinetic 
equation taking into account collective effects as in the 
3D case. This point will be further explored in the present 
paper. The kinetic theory of Coulombian interactions in 
d = 2 has been studied independently by Benedetti et al. 
[2"3"| using a different approach. They calculated the fric- 
tion force experienced by a test particle by considering 
the effect of a succession of binary encounters. The cal- 
culations are difficult because they involve the differential 
cross section which has not an explicit form in d = 2. They 
however managed to obtain asymptotic results for small 
and large velocities in the case where the distribution of 
the field particles is a step function. Our approach, based 



1 Memory effects can be important for gravitating systems 
because the temporal correlation function of the force decreases 
algebraically as t~ |17|. However, this slow decay mainly re- 
sults in logarithmic divergences of the diffusion coefficient [18) . 
so that the Landau equation can stiil be used successfufly. Fur- 
thermore, as shown in Saslaw p"], non-Markovian effects are 
important for long-range collisions but not for short-range col- 
lisions. Since long-range collisions are more gentle than short- 
range collisions, memory effects are somehow 'washed-out' in 
the complete collision process. Memory effects can be impor- 
tant for other systems with iong-range interactions, iike the 
HMF model (see below), close to the critical point T c because 
the timescale of the exponential decay of the correlation func- 
tion diverges as T -» T c |lAl2^l2ll22[ . 
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on the Landau or Lenard-Balescu equation in d = 2, does 
not require the expression of the differential cross section, 
just the Fourier transform of the potential of interaction. 
In this paper, we will calculate the diffusion coefficient and 
the friction force in the case where the distribution of the 
field particles is a step function and obtain results that 
are compatible with those of |2Sj ■ We will also discuss the 
stochastic properties of the fluctuating force created by 
a 2D Coulombian system by making a parallel with the 
approach of Chandrasekhar & von Neumann [H] for the 
gravitational force in d = 3. 

In d = 1, the Landau and Lenard-Balescu collision 
terms cancel out indicating that the collisional evolution 
of the system as a whole is due to higher order correlations 
in the 1/N expansion [22] - However, if we use this kinetic 
theory to describe the evolution of a test particle in a bath 
of field particles with any (stable) steady distribution of 
the Vlasov equation, we obtain a Fokker-Planck equation 
in which the diffusion coefficient and friction force can be 
easily calculated. This type of Fokker-Planck equations 
has been studied in |19I24I20I25| in connection with the 
HMF model. A kinetic theory of the HMF model has been 
developed in Bouchet and Bouchet & Dauxois [23 
by analyzing the stochastic process of fluctuations and 
calculating the first and second moments of the velocity 
increment (Av) and ((Av) 2 ). This is a particular case, 
for a cosine potential of interaction in one dimension, of 
the general Fokker-Planck approach developed in plasma 
physics (see Chap. 8 of Ichimaru |2S|)- These authors used 
this kinetic theory to determine the velocity correlation 
function and found an algebraic decay (see also |27I28| ) 
which is consistent with direct numerical simulations of 
the ./V-body system. A kinetic theory of the HMF model 
was developed independently by Chavanis (see discussion 
in [201) by using the projection operator formalism. This 
approach does not take into account collective effects but 
these terms can be obtained from the Lenard-Balescu the- 
ory 22 . Chavanis & Lemou used this kinetic theory 
to study the relaxation of the distribution function tails 
and showed that it has a front structure. Non-ideal effects 
in the kinetic theory (non-Markovianity, spatial inhomo- 
geneity,...) have been further discussed in 29 30 . In the 
present paper, we shall develop a kinetic theory of homo- 
geneous one-dimensional systems with an arbitrary form 
of weak long-range potential of interaction, starting di- 
rectly from the general Fokker-Planck equation given in 
|22| . which takes into account collective effects. We shall 
determine the diffusion coefficient and analyze the effects 
of a mass distribution among the particles. 

We note finally that kinetic theories of systems with 
long-range interactions have also been developed for point 
vortices in two-dimensional hydrodynamics and non-neutral 
plasmas confined by a magnetic field. Their form is related 
to, but different from, the Landau and Lenard-Balescu 
equations. The intrinsic reason is that point vortices do 
not have inertia contrary to electric charges and stars. 
Hence, the coordinates x and y are canonically conjugate. 
Kinetic equations have been derived independently by Du- 
bin & O'Neil |31l32j from the Klimontovich approach and 



by Chavanis |33j from projection operator technics. On 
the other hand, by using an analogy with stellar dynamics 
and Brownian theory, a Fokker-Planck equation describ- 
ing the stochastic evolution of a test vortex in a bath of 
field vortices at equilibrium was derived in |34I33| . The 
diffusion coefficient and the drift term are related to each 
other by an appropriate Einstein relation and they are in- 
versely proportional to the local shear created by the field 
vortices. The statistics of the velocity fluctuations arising 
from a random distribution of point vortices has been in- 
vestigated in Chavanis & Sire |35I36| by using an approach 
similar to that developed by Chandrasekhar & von Neu- 
mann [2j for the gravitational force. The analogy between 
stellar systems and 2D vortices is discussed in |57) . 

In this paper, we shall complete previous investigations 
on the kinetic theory of systems with weak long-range in- 
teractions by developing a general formalism valid for a 
large class of potentials of interaction in d dimensions and 
for multi-components systems. In Sec. 12.11 we recall the 
Lenard-Balescu and Landau equations describing the evo- 
lution of the distribution function (DF) of a system of par- 
ticles in interaction. In Sec. 12.21 we introduce the Fokker- 
Planck equation describing the evolution of a test particle 
immersed in a bath of field particles. We give the gen- 
eral expressions of the coefficients of diffusion and friction. 
In Sec. 12.31 we restrict ourselves to isotropic distribution 
functions and derive the simplified form of the Fokker- 
Planck equation. We show how its stationary solutions 
are related to the distribution of the bath and we note 
that the test particle does not in general relax towards 
the distribution of the bath except (i) if it is Maxwellian 
(thermal bath) (ii) in d = 1 for single-species particles. 
In Sec. 12.41 we consider the case where the field particles 
have a Maxwellian distribution. We find that the Fokker- 
Planck equation becomes similar to the Kramers equation 
but the diffusion is anisotropic and depends on the veloc- 
ity of the test particle. The coefficients of diffusion and 
friction are related by an Einstein relation. We derive its 
general form for a multi-components system. In Sec. 12.51 
we derive the explicit expression of the diffusion coefficient 
of a test particle in a thermal bath in d — 1, 2, 3. In Sec. 
12.61 we estimate the relaxation time of the test particle 
to the Maxwellian distribution (thermalization) and show 
that it scales as Nt d ■ We also investigate the effect of mass 
seggregation on the relaxation time. Specific applications 
are given for gravitational systems and for the HMF model 
in Sec. 12.71 In Sec. [3 we consider the case d = 3. We show 
that the results obtained by Rosenbluth et al. ^3] can be 
alternatively obtained from the Landau equation. This has 
not been noted by the previous authors and this can re- 
inforce the connection between plasma physics and stellar 
dynamics. As an application of the Rosenbluth potentials, 
we compute the diffusion coefficient and the friction force 
in the case where the distribution of the bath is a step 
function. In Sec.0] we consider the case d = 2. We obtain 
the general expressions of the diffusion coefficient and fric- 
tion force valid for any isotropic distribution of the field 
particles in terms of integrals of Bessel functions. As an 
illustration, we compute them for the Coulombian inter- 
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action when the distribution of the bath is a step function 
and we obtain results similar to those of |23| up to a fac- 
tor 2(ir — 2)/ir that may be related to the unknown large 
scale cut-off. In Sec. 14.41 we analyze the statistics of fluc- 
tuations of the Coulombian force in d — 2. We show that it 
is given by a marginal Gaussian distribution intermediate 
between normal and Levy laws. In particular, it presents 
an algebraic tail scaling as ~ F~ 4 . This is analogous to 
the distribution of the velocity created by a gas of point 
vortices in two-dimensions |35I36| . We also discuss the ori- 
gin of the linear divergence of the diffusion coefficient for 
Coulombian systems in d = 2 and how this can be cured 
by considering collective effects. Finally, in Sec. we pro- 
vide the general form of the Fokker-Planck equation in 
d = 1 and give explicit expressions of the diffusion coeffi- 
cient and friction force for any (stable) steady distribution 
of the bath and in the case where the test particle has not 
necessarily the same mass as the field particles. 



We have introduced the abbreviations /j = fi(v,t) and 
f'i = /i(v',<)- Furthermore, u = v — v' is the relative 
velocity between two particles. It is also implicitly under- 
stood that there is summation over repeated greek indices 
which denote the cartesian coordinates of the vectors. 

For systems with weak long-range interactions, the Lenard- 
Balescu equation gives the correction of order 1/N to the 
Vlasov equation which is recovered for N — ► +oo (for ho- 
mogeneous systems the Vlasov equation simply reduces to 
df /dt = 0) . The proper thermodynamic limit is defined in 
|21l22j . It is such that the amplitude of the potential of in- 
teraction scales as u ~ I /N (weak coupling) while the en- 
ergy per particle E/N , the inverse temperature (3 and the 
volume V are of the order unity. Noting that u 2 ~ 1/N 2 
and / ~ N, we find that the Lenard-Balescu collision term 
is of order 1/N, i.e. df/dt = ^Q(f) with Q(f) ~ /. In 
d = 3 and d = 2, it is easy to show that the Lenard-Balescu 
equation relaxes toward the Maxwcllian 



2 General results 

2.1 Evolution of the system as a whole: the Landau 
and Lenard-Balescu equations 

We consider a system of particles with long-range inter- 
actions whose dynamics is described by the Hamiltonian 
equations 



dr, 
dt 



dH 



dv t 
dt 



dH 
Or,, ' 



where 



H 



N 9 



i<3 



(1) 



(2) 



where u,j = it(r^ — Vj) is a binary potential of inter- 
action depending only on the absolute distance between 
particles. We assume that there exists X species of par- 
ticles {rrii, 7i}j=i x and we denote /j(r, v, t) the distri- 
bution function of particles of species i normalized such 
that J fidvdv — A^m, gives the total mass of particles of 
species i. We consider homogeneous systems and we as- 
sume that the potential of interaction is long-range and 
of weak amplitude. Then, the evolution of the distribu- 
tion function of species i is given by the Lenard-Balescu 
equation 
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dt K ' dv» 
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with the dielectric function 

e{^u J ) = l + {2K) d u{k)Y J (^ 
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Therefore, finite effects select the Maxwell distribution 
among all possible stationary solutions of the Vlasov equa- 
tion. The collisional relaxation time scales as ~ Ntn 
where tr> is a dynamical time (which can be taken of or- 
der unity). For Newtonian interactions in d = 3, there is a 
logarithmic correction hi A ~ In iV so that the collisional 
relaxation time scales as tn ~ (N/ lnN)tn- Alternatively, 
for d = 1, the Lenard-Balescu operator cancels out so that 
the collisional evolution is due to terms of higher order in 
1/AT. This implies that the collisional relaxation time of 
the system as a whole scales as tn ~ N s to with 8 > 1. 

If we neglect collective effects and take |e(k, k-v)| 2 = 1 
we obtain the Landau equation 



^ = ir{2n) d -^ I J d^'d\Lk l "k v u{k) 2 5{\L ■ u 
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3 



JJ_ 

rn ; 



mjf, 



3Q V V 



f d f'3 
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as an approximation of the Lenard-Balescu equation. An 
important remark for the following is that the expression 
||fjj of the Landau equation does not involve the differen- 
tial cross section of the interaction, but only the Fourier 
transform of the potential 2 . Furthermore, as discussed 
after Eq. (|73|l . the potential of interaction only fixes the 
timescale of relaxation so that the structure of the Landau 
equation does not depend on the potential. This is a con- 
sequence of the fact that the amplitude of the interaction 
is very small so that a weak deflection approximation is 
appropriate and yields results that are relatively indepen- 
dent on the precise form of the potential of interaction. 



The Fourier transform of the potential is just the Born 
approximation to the scattering amplitude. Therefore, the 
assumptions that are made to obtain the Landau equation 
amount to a first order perturbative approximation for the 
elastic two-body differential cross section (I thank one of the 
referees for this remark). 
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The Landau equation © can be derived from the Boltz- 
mann [2] or from the Fokker-Planck equation (HI 1 4j by us- 
ing a model of binary collisions. The Lenard-Balescu equa- 
tion @ can be obtained from the Liouville equation by 
using iterative procedures and diagrammatic methods 
It can also be obtained from the BBGKY hierarchy or 
from the Klimontovich equation |38| by using approxima- 
tions which amount to neglecting some correlations. This 
can be justified perturbatively when we consider an expan- 
sion of the equations of the problem in terms of the small 
parameter ~ 1/JV (weak coupling limit). These derivations 
remain valid for other types of potentials with weak am- 
plitude and other dimensions of space so that Eqs. J3J and 
^ have a larger domain of validity than plasma physics 
and stellar dynamics. These equations will be the starting 
point of our analysis. 



2.2 Evolution of a test particle in a bath: the 
Fokker-Planck equation 



o), (ii) 



where D^ v (j — > 0) and r? M (j — ► 0) are the diffusion and the 
friction caused by species j on the test particle denoted 0. 
With these notations, Eq. (jHJ can be rewritten 
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, dP 
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Comparing with the general expression of the Fokker- 
Planck equation 



dP _ 1 d 2 
~dt ~ 2 dv^dv" 



we find that 
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d ( (Av^ 
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= 2D^ 



(13) 
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We now consider the evolution of a test particle (m, 7) 
in a bath of field particles with prescribed distribution 
fj(v) for each species (nij,jj). It has been established 
in plasma physics (see, e.g. |2S]) that the evolution of the 
velocity distribution P(v, t) of the test particle is governed 
by a Fokker-Planck equation that takes a form similar to 
the Lenard-Balescu equation (JJJ provided that we replace 
fj(v,t) by the prescribed distribution fj(v) of the bath. 
This general Fokker-Planck equation 



(Av» 

~1aT 



(15) 



Now, using the fact that K^ v depends only on the relative 
velocity u = v — v' and using an integration by parts, we 
get 

dD» v ^ f dK^ , 



dP ,df, u(k) 2 = - V / -77- ~ m i 7,'dv' = V / K» 

dt y ' dv^ J e(k,k-v 2 ^ J 



df 



e(k,k.v)P 

dP „ df' 



' dv' v 



(1)Y1 rn^-mP^L ,(7) 
\mJ ^\mj V 3 dv y dv"J KJ 



will be our main object of interest. In the following, in or- 
der to simplify the expressions, we shall take "fi — rrii (like 
for the gravitational interaction) but the generalization of 
our results for 7$ 7^ m, is straightforward. Furthermore, 
for d = 3 and d — 2 we shall make the Landau approxima- 
tion |e(k, k ■ v)| 2 = 1 (some results that relax the Landau 
approximation are given in |22| for a thermal bath). In 
that case, the Fokker-Planck equation for the test particle 
becomes 



dP 
~dt 

where 



d 



E 



( m j f j 



'|^_ mP M) dv ' (8) 



K» v = ix{2ir) d I Wk u u(k) 2 8{\i ■ u)dk. 



(9) 



We introduce the diffusion tensor and the friction term: 
D * V = E / K^mjf'jdv' = -> °): (1°) 



(16) 



Therefore, the first moment of the velocity increment can 
be rewritten 



(Ay* 
At 



E J K^im + m^dV 



E 



m + rrij 



rf{j -> 0). (17) 



We note that because of the velocity dependence of the 
diffusion tensor, the friction terms (Av) j At and jj do not 
coincide. In particular, for equal mass particles, there is 
a factor 2 between them as noted previously On the 
other hand, if the field particles have a mass m; different 
from the mass m of the test particle, the multiplicative 
factor is (m + m/)/m. In this respect, we note that the 
frictional force (Pf r ) calculated by Kandrup (2HI with his 
linear response theory is tj, not (Av) j At. This explains 
why it differs from the calculation of Chandrasekhar [3] 
by a factor (m+mj)/m (see PP- 446). The frictional 
force 1/ is an important quantity by itself because it is the 
quantity which naturally appears in the symmetrical form 
of the Landau-Lenard-Balescu collision term where the 
diffusion coefficient is placed between the two derivatives: 
d^D^du, see Eqs. ©, ® and JEJ). 
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2.3 The isotropic Fokker-Planck equation 

If the distribution function of the field particles is isotropic, 
i.e. /j(v) = fj(v), we can write the diffusion tensor in the 
form 



1 



-Dj 



1 



-D^, (18) 



D" u = ( D l{ - 



where Dn (v) and D±(v) are the diffusion coefficients in the 
directions parallel and perpendicular to the direction of 
the test particle v. On the other hand, starting from Eqs. 
<|1U|) and and using the same type of calculation as in 
Eq. H16f) . we find that the friction term can be generally 
written 



rrij ov 



(19) 



For an isotropic distribution of the field particles, using 
Eq. I|18J1 . we find that the friction term is parallel to the 
velocity of the test particle, i.e. = rjv/v. The amplitude 
of the friction vector is given by 



* J m - 



dD \\ 

dv 



d - 1 f D _Jh 



U - o). 

(20) 



If the velocity distribution of the test particle is itself 
isotropic, i.e. P(v, t) — P(v, t), we can rewrite the Fokker- 
Planck equation l|12fl in the form 



dP 
~dt 



1 d 

,jd-l Q v 



(21) 



where we have used 



D^v" = (D\\ - 



d-1 



d-1 



fli / = D\\v>- 



Equation (|21|l can also be written as 
1 d 



dP 
~dt 



„d-l 



dv 



, , i dP dU 



where we have introduced the effective potential 



U(v) 



D\\(V) 



-dv'. 



(22) 



(23) 



(24) 



Equation (|23[1 relaxes towards a stationary solution of the 
form 



P eq {v)=A e - u(v \ 



(25) 



provided that this distribution is normalizable. The evo- 
lution of the high velocity tail of the distribution function 
that is solution of a Fokker-Planck equation of the general 
form 1)23(1 has been studied in |25| . 



Consider single species systems. For d — 2,3, the only 
stationary solution of the Lenard-Balescu equation © is 
the Maxwellian J5| ■ This implies that the stationary solu- 
tion of the Fokker-Planck equation Q will be equal to the 
distribution of the bath /(v) only if this distribution is the 
Maxwellian. Otherwise, P eq (v) is not equal to the distri- 
bution /(v) of the bath. We shall give an explicit example 
in Sec. 13.31 By contrast, for d = 1, the Lenard-Balescu op- 
erator cancels out for any distribution. This implies that 
the stationary solution of the Fokker-Planck equation (0 
is always equal to the distribution of the bath /(v), even 
if it is not the Maxwellian (see SecEJ). 



2.4 The Einstein relation for isothermal systems 

If the velocity distribution of the field particles is the 
Maxwellian (thermal bath) 



then 



fr- 
ofi 

dv' v 



(26) 



(27) 



Substituting this relation in Eq. (|ll|l and using the rela- 
tion 



K^v u = - u v ) = K^v", 



(28) 



which results from the identity K^ u u v — (K^ u is the 
projector in the direction perpendicular to u), we get 



if = fimD^v". 



(29) 



This can be viewed as the general expression of the Ein- 
stein relation in our context (this is the most general re- 
lation that Eq. (|12fl must satisfy in order to admit the 
Maxwell distribution as a stationary state). In that case, 
the Fokker-Planck equation l(12J) can be written in a form 
similar to the Kramers equation |41| , but with an anisotropic 
tensor depending on the velocity of the test particle: 



dP 
~dt 



d 

dvf* 



( BP 
\dv v 



[3mPv L ' 



(30) 



The stationary solution of this equation is the Maxwellian 

P eq {v) ^Ae-P" 1 ^. (31) 

Note also that according to Eq. I)22|) we have 

q = /3mZ>||V. (32) 

Therefore, if we assume that the velocity distribution of 
the test particle is isotropic we get 



dP 



1 d 

jd-l Q v 



dP 



(33) 
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If we momentarily neglect the velocity dependence of the 
diffusion coefficient, this Fokker-Planck equation can be 
obtained from a Langevin stochastic process of the form 



dr dv 



2£>||R(i), 



(34) 



where R(t) is a white noise satisfying (R(f)) = and 
(Ri(t)Rj(t')) = 5ijS(t — t') and the friction coefficient is 
given by the Einstein relation £ = 0mDu. This relation is 
necessary to obtain the Maxwellian distribution at equi- 
librium and this is why it is rather independent on the 
details of the microscopic model. 

As indicated previously, since the diffusion coefficient 
depends on the velocity, r) is not exactly the friction force 
so that Eq. I|32() is not the proper form of Einstein relation. 
Using Eq. I|17f) we have 



At 



(35) 




Fig. 1. Diffusion coefficients Gu(x), G±(x) and friction force 
xGu(x) for a thermal bath in d — 3. 



If the field particles all have the same mass m f we obtain This can be rewritten 



At 



-f3(m + m/)£)iiv, 



(36) 



so that the proper friction coefficient is 

? = P(m + mf)D t (37) 

This is the proper form of the Einstein relation in the 
present context. Note that it involves the sum of the mass 
of the test particle and of the field particles. Using 



- 2 m f {v ) f = -k B T=- ) 
the Einstein relation 1)37(1 can be rewritten 



((^||) 2 ) _ 2 m f 



£' At dm + rrif 



(v 2 )f- 



(38) 



(39) 



2.5 The diffusion coefficient for isothermal systems 

Inserting the identity 



S(x) 



. dt 
2V 



(40) 



in Eq. (|5J) and performing the Fourier transform on v' of 
the Gaussian distribution l|26|) in Eq. l]lUp. we find after 
another Fourier transform on t that the diffusion tensor 
can be expressed as 



= 7r(27r) 







2tt 



1/2 



3/2 



/ dk'^-j—uikfe- 13 " 1 !^ 



(41) 



a \ V 2 

P \ sr^ _ _ 3/2 



2;t,/ 

, / k d u{kfdkG^ ' 



. lv ( 0m j 



where 



(42) 



(43) 



where we have noted k = k/fc. In d = 3, introducing 
a spherical system of coordinates with the z-axis in the 
direction of x, we obtain 



1 x^x u 1 
= (G„ - ^G ± )— — + -G ± 5^ 



(44) 



with 



2tt 3 / 2 2vr 3 / 2 
G|, = G{x), G± = [*(a;) - G(x)}, (45) 



where 



G{x) = ^ L te dt= — ^ — • 



and <P(x) is the error function 

9 r x 



(46) 



(47) 



In d — 2, using a polar system of coordinates with the 
cc-axis in the direction of x, we get 



G"" = (G|, - G ± )^4- + 



(48) 
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d=2 










Gpe,p(><) 




7 \V XG pB.a( X ) 







8 10 



Fig. 2. Diffusion coefficients G||(a;), G±(x) and friction force 
xG\\ (x) for a thermal bath in d — 2. 



1.5 - 



1 - 



0.5 - 




Fig. 3. Diffusion coefficient G(x) and friction force xG(x) for 
a thermal bath in d — 1. 



with 



Gii = 7re 2 



G i = 7re 2 



Jo 



h 



2 \ 1 



Finally, in d = 1, we obtain 
G(x) = 2e' 



(49) 
(50) 

(51) 



The diffusion coefficients and friction force are plotted in 
Figs. 11131 for different dimensions of space. 



2.6 The relaxation time 

We can use the above results to estimate the relaxation 
time of the test particle towards the Maxwellian distribu- 



tion (thermalization). We consider the relaxation of a test 
particle with mass m in a thermal bath of field particles 
with mass mj. If we set x = y/ (3m f /2v, the Fokker-Planck 
equation 1)30(1 can be rewritten 



dP 
~dt 



1 d 



G» v (x) 



dP 



t R dx^ 

where tR is a reference time given by 



2^Px» 
m f 



(52) 



1 



V2 . d pm f 



d 3 / 2 (2ir) 



rnf 



+ DC 



k d u(k) 2 dk, (53) 



where = dj (/3m f) is the mean-squared velocity of the 
field particles. Introducing the dimcnsionless function 

T)(k) = -(2Tr) d u(k)(3m f p, (54) 

we can rewrite the reference time (|53|) in the form 



1 



\8d) 



1/2 



K d r){n/L) 2 dK, (55) 



where L — V x l d is the size of the system and k — kL is 
dimcnsionless (n — p/nrif is the numerical density of the 
field particles). We define a typical dynamical time by 



to = 



Vmf 



(56) 



Then, introducing the number N = nL d of field particles, 
we can finally write the reference time in the form 



tR — CdNtn, 



(57) 



where 



1/2 1 



(27T) d 



+ OC 



K d r](K/L) 2 dn, (58) 



is a dimensionless number. 

We can get an estimate of the relaxation time by the 
following argument. If the diffusion coefficient were con- 
stant, the typical velocity of the test particle (in one spa- 
tial direction) would increase like 



>vf> 



2D\d. 



(59) 



The relaxation time t r is the typical time at which the 
typical velocity of the test particle has reached its equi- 
librium value (v 2 )(+oo) = d/(m(3) = (jnjlm^v 2 ^^ so that 

((Av) 2 )(t r ) = (v 2 )(+oo). Since D\\ depends on v, the de- 
scription of the diffusion is more complex. However, the 
formula 



1 mt 



d m 2D\\(v m f) ' 



(60) 



resulting from the previous arguments with Du = D» (v m f) 
should provide a good estimate of the relaxation time. Us- 
ing Eq. 142(1 and comparing with Eq. ((55(1 we obtain 



mt 

Kd tR, 

m 



(61) 
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where K d = l/[AG\\{^/d/2)). We find K 3 = 0.13587547..., 
K 2 = 0.16286327... and K x = 0.20609016.... We can also 
estimate the relaxation time by t' r — where f is the 
friction coefficient. Using the Einstein relation £ = Duftm 
with _D|| = D\\(v m f) we find that 



t r — 2i r . 



(62) 



Combining the previous results we find that the relax- 
ation time scales as 



N — -to- 
rn 



(63) 



We expect that this result also yields a good estimate of 
the relaxation time of a test particle towards the station- 
ary distribution (|25|l in the case of a non-thermal bath. In 
that case, we must justify that the DF of the bath does 
not change on that interval (see below) . We come therefore 
to the following conclusions. For equal mass particles, the 
relaxation time of a test particle in a bath scales as Nto- 
For d = 2,3, the relaxation time of the system as a whole 
also scales as Nto (see Sec. l2.1|) . Therefore, a non-thermal 
bath will change on this timescale. Only the Maxwellian 
distribution is stationary on the timescale Nt u ■ Thus, for 
m = mf, the test particle approach can be developed 
only for a thermal bath. Now, consider the relaxation of 
a test particle with mass m in a bath of field particles 
with mass to/. In that case, the relaxation of the test par- 
ticle is changed by a factor m//m. If to 3> to/ we find 
that the relaxation of the test particle (~ (m/ /m)Ntn) 
towards Eq. (|25|l is much faster than the relaxation of the 
system of field particles as a whole (~ Ntn) towards the 
Maxwellian f eq . Therefore, in that limit, when we focus 
on the evolution of a test particle, it is possible to con- 
sider that the distribution of the field particles /(v, t) is 
"frozen" even if it does not correspond to statistical equi- 
librium. This is because /(v, t) evolves much slower than 
P(v,t). This remark justifies to consider equations of the 
form 10 with / ^ f eq . In astrophysics, the case m^> m,f 
could be relevant to describe the stochastic evolution of a 
black hole at the center of a galaxy or the dynamics of a 
globular cluster (m ~ 10 6 to/) passing through a galaxy. 
Note that in d = 1 the relaxation of a test particle in 
a bath is given by Eq. (|63l) while the relaxation of the 
field particles as a whole towards statistical equilibrium is 
larger than Ntn since the Lenard-Balescu collision term 
cancels out in d = 1. Therefore, in that case it is justified 
to consider the relaxation of a test particle in a bath with 
any distribution f(v) (stable with respect to the Vlasov 
equation) even if to = to/ (see Sec. EJl- 



2.7 Examples 

Let us apply these results to particular systems considered 
in |21I22| . For the gravitational interaction in d = 3, the 
Fourier transform of the potential is 



(27r) 3 u(fc) 



4ttG 



(64) 



The reference time l)53|l can be written 



t R = 0.482 



J mf 



nm 2 G 2 In A ' 



(65) 



where lnyl is the Coulombian factor |15|. Using Eqs. (|54|1 
and ((Mil , we get 



(66) 



where kj = {AirGflpmf) 1 / 2 is the Jeans wavenumber. Then, 
Eqs. |57J)-lEBl) lead to 



11.3 N 
r\ 2 In A 



(67) 



where we have defined 77 = /SGNmj /R with L 3 = V 
(4/3)tt i? 3 and we recall that ln^L - IniV. 
For the HMF model in d = 1, we have 

k 

U n = —-7-(S n ,l + S n -l). 

The reference time (1531 can be written 



(68) 



t R = 40.1 



nm 2 k 2 ' 



(69) 



where we have replaced / +O ° dk by X^n^o m ^ ne discrete 
case. Using Eqs. (fS^ and (|rJ5|) . we get 

f3Nm 2 f k 

Vn = V(S n ,l + S„ ,-1), = (70) 



4tt 



Then, we obtain 



tR 



0.254 



-Nt D . 



(71) 



For the HMF model, the relaxation towards the Maxwellian 
is not exponential due to the rapid decay of the diffusion 
coefficient |24l25j so that this reference time only gives a 
timescale of relaxation. 



3 The case d = 3 

3.1 Rosenbluth potentials 

In d = 3, it is possible to obtain a simple expression of the 
coefficients of diffusion and friction, expressed in terms 
of elementary integrals, for any isotropic distribution of 
the bath (not only the Maxwellian). To that purpose, it 
is useful to introduce the so-called Rosenbluth potentials 
|14ll5l39j . Introducing a spherical system of coordinates 
with the z-axis in the direction of u, we find after elemen- 
tary calculations that Eq. © can be written 



A. 



5^u 



(72) 
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where 



A, = 8tt 5 



k 3 u n ^ 2 



(kydk. 



(73) 



In particular, for the gravitational potential we find that 
A3 = 2nG 2 \nA where \nA — J + °° dk/k is the Coulom- 
bian factor which has to be regularized with appropriate 
cut-offs (see the Introduction). Using Eq. (|72|l . the Landau 
equation © can be written 



dt 



A- 



d_ 



S^u 2 



u 3 \ dv u 



dv' 1 



(74) 



with A = 2irG 2 mlnA. Equation l|74|) is the original form 
of the Landau equation (we have used the notations of as- 
trophysics but the connection with the notations of plasma 
physics is obtained by making the substitution Gm 2 «-* 
— e 2 ). It applies to weakly coupled systems. We note that 
the potential of interaction appears only in the constant A 
which merely determines the timescale of relaxation. The 
structure of the Landau equation is independent on the 
potential. Using the identity 



d 2 u _ S^u 2 - u^u v 
dv^dv v u 3 

we can write the diffusion tensor 110H in the form 



D> 



As. > m 



d 2 9l 



^ " ~ 3 dv»dv v 
3 



(v), 



where 



ff(v)= / /(v')|v-v'|dv'. 



(75) 



(76) 



(77) 



On the other hand, integrating by parts and using the fact 
that K^ v depends only on the relative velocity u = v — v' , 
the friction term 111(1 can be written 



= — m 



Using the identity 



-2 A,— = 2 A, 



d (\ 



we can rewrite the friction term in the form 



where 



(78) 



(79) 



(80) 



(81) 



The functions <?(v) and /i(v) are called the Rosenbluth 
potentials. Using Eqs. (|7f>|) and l(8()|l . we can readily check 



that the first and second moments of the velocity incre- 
ments given by Eqs. (|14|) and (|17|) return the expressions 
obtained by Rosenbluth et al. 14 and that the kinetic 
equation derived in ^1] is equivalent to Eq. 1)74(1 , although 
written in a different form. Therefore, we have recovered 
the kinetic equation of stellar dynamics directly from the 
Landau equation. This shows that the Landau equation 
is equivalent to the kinetic equation derived by Chan- 
drasekhar |3] and Rosenbluth et al. starting directly 
from the Fokker-Planck equation and evaluating the first 
and second moments of the velocity increments by consid- 
ering a succession of binary encounters. This was expected 
because the Landau equation carries the same type of as- 
sumptions. However, it is surprising that the relation to 
the Landau equation was not mentioned in [3114115) . In 
particular, these authors write the Fokker-Planck equa- 
tion in the form ((13(1 with the diffusion coefficient in the 
second derivative (d ll d v D 111 ') while a more symmetric form 
is the Landau equation ((74(1 where the diffusion coeffi- 
cient is inserted between the first derivatives (d^D^ d u ), 
see Eq. ((12(1 . It is interesting to note, for historical rea- 
sons, that this symmetric form (from which we immedi- 
ately deduce all the conservation laws of the system and 
the i?-theorem 0]) has escaped to the study of stellar dy- 
namicists [3114115) while the Landau equation was known 
long before in plasma physics. 



3.2 Diffusion and friction 

We shall see that the Rosenbluth potentials can be eas- 
ily calculated in d = 3 when the field particles have an 
isotropic velocity distribution. The results of the previous 
section remain valid in d = 2 provided that we replace A3 
by A2 defined in Eq. I(122() . However, the calculation of 
the Rosenbluth potentials is apparently more complicated 
in d = 2 than in d = 3 and this is why we shall use an- 
other method in d = 2 to obtain the diffusion and friction 
coefficients (see Sec. 14.11) . 

If the field particles have an isotropic velocity distri- 
bution, the Rosenbluth potentials in d — 3 take the par- 
ticularly simple form [14115139] : 



h(v) 



An 



f{vi)vldvi 



fiv^vidvi 



,(82) 



. . Attv 
9(v) = — 



r+00 



3v 3 \ 
— - + vvi f{v\)dv\ 

V J 



(83) 



Furthermore, when g = g(v) the diffusion tensor ((76(1 can 
be put in the form of Eq. ((18(1 with 



9j 



dv 2 



Di =2A : 



£ — ' J 11 rlii 



(84) 



(85) 
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Using Eq. H83[l we obtain 



+v 



vifj(v 1 )dv 1 



(86) 



D 



-2v 



vifj(vi)dvi 



(87) 



On the other hand, when h = h(v) the friction term l|8U|l 
can be written 




-2A 3 m^2 



1 dhj 

v dv 



Using Eq. we obtain 

v f v 
i/=8vrA 3 TO— V/ fj(vi)vfdvi. 

V 4 JO 



(88) 



(89) 



Fig. 4. Diffusion coefficients D«(v)/ D*, D±(v) / D* and fric- 
tion force rj(v)/ri* as a function of v/vo for a water-bag dis- 
tribution of the bath in d = 3. The normalization factors are 
D t = ^-A 3 mfr]oVo and n* = ^-A 3 mr} v . 



We note that these expressions are valid for any isotropic 
distribution of the held particles. If we substitute Eqs. I|86[) 
and (|89|) into Eq. I|21ll , we get the Fokker-Planck equation 
describing the evolution of a test particle in a bath with 
prescribed distribution fj(v). Alternatively, if we come 
back to the original Landau kinetic equation 10, assume 
an isotropic velocity distribution and substitute the gen- 
eral expressions (|8T?|l and (f%9"|) with now fj — fj{v,t) we 
obtain the integro-differential equation 



and for v > vq: 



8tt „ 
V = -T^-Asmriov. 



8tt 



D n = —A 



15 v a 



D 



A, 



5v< 



- 8ttA 3 



dt 



v 2 dv 



mj dfj / 1 
3 dv \ v ./n 



vff J (v 1 ,t)dv 1 



^ A 4 
3 v z 



(93) 

(94) 
(95) 
(96) 



vifj(vt,t)dvi +iriifi / fj(v!,t)vldvi 



(90) 



These quantities are plotted in Fig. 0] We note that the 
friction term and the diffusion coefficient in the direction 
parallel to the direction of the test particle are related to 
each other by 



describing the evolution of the system as a whole. 



2mD 



V 



m f (v 2 



h 2 ) 



3.3 Water-bag distribution 

As an illustration of the previous formalism, we shall com- 
pute the diffusion coefficient and friction force of a test 
particle with mass m when the distribution function of 
the bath (composed of particles with mass m/) is a step 
function: /(v) = i] Q for v < v and /(v) = for v > v 
with ?/o = 3p/(47twq) (p is the mass density of the field par- 
ticles). Using Eqs. tfHgft . l|87)l and (HH|) we get for v < v : 



5mD 



T v > 



V, (V < V ), 



(v > Vq). 



(97) 



(98) 



These expresssions can be compared to the Einstein re- 
lation 132fl . We note that the role of the temperature is 
played here by rrifv 2 . 

Substituting these relations in Eq. I|24[l we get a Fokker- 
Planck equation (|23[) with an effective potential 



47T 

^11 = T A 3 m fno 



(91) 



U(v) 



o m 
3mr 



-In 



3^ 
2vl 



(v<v ), (99) 



87r A ( 2 « 

Dj_ = Y A s m fVa I v - — 



(92) 



5 m v 2 

U{V) = - [ v > v 0)- 

2m f v 2 



(100) 
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The stationary solution of the Fokker- Planck equation is 

{v<v ) (101) 



P eq (v) = Ae 



{v > v Q ) 



(102) 



where A is a normalization constant. We explicitly check 
on this example that the stationary velocity distribution 
of the test particle is different from that of the bath. As ex- 
plained at the end of Sec. 12. 61 this is just a quasi-stationary 
distribution valid on a timescale (mj /m)Ntn -C Nto for 
m^> rrif because on longer timescales (~ Ntn) the bath 
distribution will change. 

We can also compute the diffusion coefficient and the 
friction force from Eqs. JBHJl, lf%7|) and lf%9"|) for a Maxwellian 
velocity distribution of the field particles. In that case, we 
directly check that the Einstein relation (|32|) holds and 
that the diffusion coefficients are given by Eqs. i|45J) and 
H4(j|) . For comparison with the water-bag model, we give 
below the asymptotic expressions of the diffusion coeffi- 
cient and friction force for a thermal bath. For v — > +oo, 
we have 



Du = 2 A : 



D± = 2A a pm f -, 

v 



rj = 2A 3 pm—. 



(103) 



(104) 



(105) 



For v — > 0, we get 



Du 



-A, 



P\ 1/2 3/2 A 3 a 
P m f [ l ^T^P rn f v 



2?r 



10' 



(106) 



D 



3 \2ir 



1/2 



3/2 



1 

10' 



1 — —/3mfV 2 



"=3* 27 



1/2 



3/2 „ 

pmy pmv. 



(107) 



(108) 



4 The case d = 2 
4.1 General expressions 

We shall provide here general expressions of the diffusion 
coefficient and friction force in d = 2 for any isotropic 
velocity distribution of the field particles. We shall use a 
method different from that exposed in Sec. |3 for d = 3. 
The following method also works in d = 3 as shown in 
Appendix [X] 

Using Eq. JSJ , the diffusion coefficient i|lU|) can be writ- 
ten 

= Tr(2TT) d J2 m 3 J dkdv'k^k'uikfSik ■ u)fj(v'). 



(109) 



Inserting the identity (|40|l in Eq. I|109|) . we get 



D» v = (2vr) 2d J2 m J J dt J dkfc^fc^fc) V k v 7j(kt), 



where we have introduced the Fourier transform 

dx 



/(k) = /(x)e 



-ik-x 



(2ir) d ' 



(110) 



(111) 



For an isotropic velocity distribution, /(k) = <P(k) de- 
pends only on the modulus of k (see the explicit expression 
below). Setting t — kt we finally obtain 



D"" = (2n) 2d J2 



n, l I dr / dk- 



-u(kye tk - VT ^(T) 
(112) 



where k = k/fc. This expression is valid in d dimensions. 
We now specialize on the case d = 2. Introducing polar 
coordinates with the x-axis in the direction of v and using 
the identities 



ix cos 9 



dQ = 2irJ (x), 



e lxcos9 sin 2 ed9 = 2n^-. 



(113) 
(114) 



we find that D^ v is given by Eq. l |T8|l with 

1 f +OQ 

Dp = (2tt) 5 Yl m J- J k 2 u{k) 2 dk 



f 

Jo 



dx 



J (x) 



Ji{x) 



(115) 



Interestingly, these are the same asymptotic behaviors (same 
scaling) as in the water-bag model. For recent studies con- 
cerning the motion of a "test particle" in a thermal bath 
in relation with space plasmas, and for some numerical so- 
lutions of the corresponding Fokker-Planck equation see, 
e.g., Shizgal 03| an d references therein. 



1 f H 

D± = (2ir) 5 m i ~ J k 2 u{kfdk 

x \vJ 



(116) 
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Now, in d = 2, 



D 



2nA 2 m f r ]0 v R 1 _ (^) , (126) 



1 



#(*0 = 7T / fMMkv^dvL (117) 
27r Jo 

Substituting this expression in Eqs. (|1 15|) and l|116[) and 
introducing the notations 



J (x) 



Ji(x) 



J (Xx)dx, (118) 



+ OC 



Ji(x) 



Jo(Xx)dx, 



(119) 



we obtain 



D|, =2ii 2 ^m ] l| + (^fidwi, (120) 



27rA2 m, — 



/ + /jW-flf-)fi*i,(121) 



where 



Sir* 



+ 00 



kru{kydk. 



(122) 



Finally, the friction is given by Eq. I|2U|I . Note that these 
expressions are valid for an arbitrary isotropic velocity dis- 
tribution of the field particles. They remain true if the dis- 
tribution function depends on time. Therefore, by substi- 
tuting these results in the Landau equation © we obtain 
a self-consistent kinetic equation for an isotropic distri- 
bution function f(v,t) which is the counterpart in d — 2 
of Eq. I|9(J|I in d — 3. However, because it involves com- 
plicated functions (|118fl and I119f) . its expression is less 
explicit. 



4.2 Water-bag distribution 

The previous relations are quite general. Let us now con- 
sider for illustration the case where the distribution func- 
tion of the field particles (with mass m/) is a step func- 
tion: f(v) — r/o for v < vo and f(v) = for v > vq where 
7/0 = p/(ttVq). Using the identity 

xJ (x) = [xJx{x)]' : (123) 

the Fourier transform l|117|) of the distribution function is 



2tt£ 



Jidvo). 



(124) 



Inserting this expression in Eqs. (|115fl and i|116|) . we ob- 
tain 



Dn 



2TrA 2 mfri v Q R l \ ) , 



(125) 



where we have introduced the functions 



iZ.,(A) 



dx 



J Q (x) 



Ji(x) 



-Ji(Ax), (127) 



R±W 



°° dx il^lj^Xx). 



to x 
Using Eq. (|20|l . the friction term is given by 



(128) 



V 



m 



dD 



dv v 



(129) 



The integrals (|127fl and 1)1281) can be expressed in terms 
of hypergeometric functions. Their asymptotic behaviors 
are derived in Appendix 151 Using these results, we obtain 
the asymptotic behaviors of the diffusion coefficients and 
friction force. For v — > +oo, we get 



1 , 4 1 

D\i = -irA 2 m f r] v Q -^ 



D± = nA 2 mfr]ov -, 



nA 2 mri v ^T. 



We note that the leading term 



(130) 



(131) 



(132) 



of the friction 77 is 



due to the diffusion in the direction perpendicular to the 
velocity of the test particle, i.e. the term D± in Eq. i|129|) . 
The contribution of the diffusion coefficient in the direc- 
tion parallel to the velocity of the test particle decreases 
more rapidly like u -4 . For v — > 0, we get 



Dn 



TrA 2 mfr] v 



8 \v 



Dj_ = irA 



2 mf7j v 



1 - 



8 \v 



77 = irmA 2 r] — . 



(133) 



(134) 



(135) 



On the other hand, using Eqs. I|129|) and l|219|) . we find 
that the friction coefficient diverges for v — > vq like 

7/ = 2A 2 mf] (2 - 27 - ln2 - 2-0(3/2) — In 1 1 - v/v \). 

(136) 

The diffusion coefficients and friction force for a water-bag 
distribution of the bath are plotted in Fig. |3J 
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Fig. 5. Diffusion coefficients Du (v)/D*, D±(v)/ D* and fric- 
tion force rj(v)/rj, as a function of v/vo for a water bag dis- 
tribution of the bath in d = 2. The normalization factors are 
D* = 2nA2mfrjoVo and ry» = (m/mf)(D*/vo). 



We can also use the previous formalism in the case 
where the distribution of the bath is Maxwellian (thermal 
bath), in which case 



*i(0 



Pi 



(27T) 



2,-!n 



(137) 



Substituting this expression in Eqs. fTT5|l and fTTnj) and 

carrying out the integrations, we can show that they re- 
turn the results l|4Tj|l and (|5*tj|l obtained by a different 
method. For comparison with the water-bag model, we 
give below the asymptotic expressions of the diffusion coef- 
ficient and friction force for a thermal bath. For v — > +oo, 
we have 



D« = An 



1 

'■0v* 



D± = A 2 prrif 



(138) 



(139) 



V 



1/2 



3/2 fl 

pmy pmv. 



(143) 



Interestingly, these are the same asymptotic behaviors (same 
scaling) as in the water-bag model. This is also the case 
in d — 3 (see Sec.|3J). 



4.3 Example: Coulombian plasma 

In a recent paper, Benedetti et al. [23j have considered a 
2D model of Coulomb oscillators interacting via a poten- 
tial of the form 



1 

N' 



In 



(144) 



They have calculated the friction force K = ((Av))/At 
experienced by a test particle by evaluating the average 
variation of the velocity caused by a succession of bi- 
nary encounters. The calculations are relatively difficult 
because the explicit expression of the differential cross sec- 
tion for a Coulombian interaction is not known in d = 2. 
We shall reconsider this problem with our approach based 
on the Landau or Lenard-Balescu equation which does 
not require the explicit expression of the differential cross 
section, but only the Fourier transform of the potential 
of interaction. As noted after Eq. (|74J) . for weakly inter- 
acting systems, the form of potential just determines the 
timescale of relaxation through the constant A. 

First, we note that in d = 2 the potential (|144|l is 
solution of the differential equation 



Au=- 2 -f5{v). 



This immediately leads to 



(27r) 2 u(fc) 



2?r£ 



(145) 



(146) 



Using Eq. (|146|) . we find that the constant l|122[) is given 

by 



a 2 < 2 , 
A 2 = — irA 

2 N 2 



(147) 



For v — > 0, we get 



7T , / 
Du = ~A 2 — 
11 2 V 2tt 



D, = -A 2 



i] = A 2 pm^r. 



1/2 



3/2 A 3 2 

pnif \l——pmfV + 



1/2 



3/2 

prrif 



-0mfv 2 



(140) 



(141) 



(142) 



where 



50 dk 



A= 



(148) 



We note that for a Coulombian potential in d = 2, the 
Coulomb factor A diverges linearly with the distance for 
A = 27r/fc — > +oo. This contrasts from the 3D case where 
the divergence is only logarithmic. In paper I22j . we have 
suggested that this divergence would be cured (as in the 
3D case) by using the complete form of the Lenard-Balescu 
equation including the collective effects encapsulated in 
the dielectric function. These calculations are completed 
in Appendix El However, in a first approach, we shall ne- 
glect collective effects and argue phenomenologically that 
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Fig. 6. Friction force as a function of the velocity for an 
isothermal distribution of the bath (theoretical curve obtained 
with our approach). For comparison, we have adopted the same 
normalization as in |2o) . 



Fig. 7. Friction force as a function of the velocity when the 
distribution of the bath is a step function (theoretical curve 
obtained with our approach). For comparison, we have adopted 
the same normalization as in |23| . For N —* +oo, the friction 
force diverges logarithmically at the velocity vq. 



A should scale like the Debye length in 2D. Note that there 
can be a numerical factor between A and the Debye length 
so that our approach will only provide an estimate of the 
diffusion coefficient. 

If we specialize on the case where the distribution of 
the bath is a step function, we can use the results of Sec. 
IQwith 



N 



(2^) 2 k B TR 2 ' 



2^/kgT, (149) 



where, following |23|. we have introduced a typical ra- 
dius via p — N/(ttR 2 ) and defined the "temperature" via 
(v 2 ) = 2k B T (we have taken m = 1). Recalling that for 
identical particles the friction force is given by K = 2r\ 
and using the asymptotic results ((132(1 and 1(135(1 . we get 
for v — > +00: 



K 



e 



NR 2 



A 



1 



and for v — > 0: 



K = 



1 e 



2 NR 2 (k B T) 3 / 2 



(150) 



(151) 



These asymptotic results agree with those obtained by 
Benedetti et al. [231 with a different method, up to a factor 
2(n — 2)/n ~ 0.727.... Owing to the remark following Eq. 
1)148(1 . we should not give too much credit on the exact 
value of the numerical constant (see Appendix [0 for a 
more precise determination of the value of A). However, it 
would be important in future works to determine whether 
our approach and that of 23 are really equivalent or not. 

We emphasize that our approach provides the expres- 
sion of the diffusion tensor D^ v and friction force rf for 
any value of the velocity; these quantities have not been 



calculated in . When the distribution of the field par- 
ticles is Maxwellian (thermal bath), they are given by the 
analytical formulae (|49|l and (|50|l obtained in 22|. The 
friction force obtained from Eq. 1(32(1 is plotted as a func- 
tion of the velocity in Fig. We have multiplied our ex- 
pression by the factor 2 (71" — 2)/ir. In that case, our analyt- 
ical formula of K matches remarkably well, in the whole 
range of velocities, the curve obtained numerically by |23| . 
For a water-bag distribution of the bath, the diffusion co- 
efficients and friction force can be expressed in terms of 
integrals of Bessel functions ((125(1 - 1(129(1 or, equivalently, in 
terms of Hypergeometric functions using Eqs. ((211(1 . 121211 
and 1(218(1 . The friction force obtained from Eq. 1(129(1 is 
plotted as a function of the velocity in Fig. [7| Again, we 
have multiplied our expression by the factor 2(7r— 2) / 7r and 
we get an excellent agreement with the numerical results 
of [231 - We note that in our approach valid for N — > +00, 
the friction force diverges for v — > vq. Using Eq. 1(136(1 . we 
find that the divergence is like 



K. 



e 



A 



NR 2 k B T 



6 In 



2^/k^T 



(152) 



where a = (2 - 2 7 - In 2 - 2^(3/2))/tt ~ 0.025287 and 
b = 1/tt ~ 0.3183. This behavior is consistent with the 
results of |23| who obtain the K — v curve from A^-body 
simulations and find that the maximum velocity increases 
as N increases. For finite N systems, there is no singu- 
larity. The singularity appears for N — > +00. Formally, 
the curve of Fig. [3 and the divergence of the friction force 
at v = vq when N — > +00 share some analogies with 
the divergence of the specific heats in second order phase 
transitions (e.g., the curve of Fig. [7| looks similar to the 
A-transition in superfluid Helium). 
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4.4 Statistics of fluctuations 



which behaves like 



In order to better understand the previous results, and in 
particular the linear divergence of the diffusion coefficient, 
we can analyze the statistics of fluctuations of the force 
produced by a random distribution of charges in d = 2. 
The force created by the charges at some point of origin 
is given by 



N 

2—1 1 



(153) 



We shall assume that the charges are randomly distributed 
with uniform distribution. In that case, the force F fluc- 
tuates and we have to determine its distribution W(F). 
Formally, this mathematical problem is the same as the 
one considered by Chavanis & Sire in their investiga- 
tion of the statistics of the fluctuations of velocity V cre- 
ated by a random distribution of point vortices in d = 2. 
Therefore, we can apply their results just by making the 
substitution F «-> V (we also note that to leading order 
(F 2 ) ~ N x (1/N 2 ) ~ 1/N so that the proper scaled vari- 
able is x = FN 1 ' 2 ). To have exactly the same notations, 
we set £/N = 7/(271"). It is shown in that the distri- 
bution of F is a marginal Gaussian distribution 



W(F) 



ri7 2 In N 



(F -C F crit (N)), (154) 



W(F) 

F crit (N) ~ 



717 



(F ^> F crlt (N)), (155) 
\ 1/2 

^-L\nN) ln 1/2 (lniV), (156) 

47T / 



4tt 2 F 4 
2 



where n is the spatial distribution of the particles assumed 
to be homogeneous. This distribution is intermediate be- 
tween Normal and Levy laws: the core of the distribu- 
tion is Gaussian but the tail decreases algebraically as 
F~ 4 . It is dominated by the contribution of the nearest 
neighbor . Note that the variance of the scaled variable 
X = FN 1 / 2 diverges logarithmically with N since 



(F 2 



717 

47T 



In TV. 



(157) 



Analytical results for the spatial correlations of the force 
are derived in |36| . Here, we shall concentrate on the basics 
in order to make the connection with the kinetic theory 
developed previously. If we neglect collective effects, the 
spatial correlation function of the force can be written (see 
Eq. (100) of Ref. |ZH): 



(F(0) • F(r)) = (27r) 2 n / k 2 u(k) 2 e 



\2 -ik-: 



dk. 



(158) 



Introducing a polar system of coordinates with the x axis 
in the direction of r, using (27r) 2 u(fc) = 7/fc 2 and intro- 
ducing a large scale cut-off A, we obtain 



(F(0) 



^ ^ 27T j i/A 



■h(kr) 



dk, 



(159) 



<F(0).F(r))^^lnf- 



(160) 



for r/A — » 0. An alternative derivation of this result is 
given in Ref. |36) . Appendix D. In particular, we find that 
the correlation function diverges logarithmically for r — > 0, 
in agreement with the result l|157[) [note that in d = 2 
the length scales as L ~ (N/n) 1 / 2 which yields the factor 
i In TV] . We shall now see how collective effects can remove 
this logarithmic divergence. The expression l|158fl neglects 
the correlations between particles. Now, applying Eq. (51) 
of Ref. [21] to the present context, it is found that the 
spatial correlation function is solution of the differential 
equation 



Ah-k z D h = /37<5(x), 



(161) 



where fc/j = (/Jrry) 1 / 2 is the Debye wavenumber. The spa- 
tial correlation is then given by 



K x ) = -7T-K a (k D x), 

Z7T 



(162) 



which is the Debye-Hiickel result in 2D. The Fourier trans- 
form of the correlation function can be written 



(2ir) 2 nh(k) 



7.2 
h D 



k 2 



fcf. 



(163) 



and the correlational energy (see Eq. (57) of paper |21j) 
can be written W c = -wyV/(4>rr)['y E +(l/2) ln((3n-f)-\n2] 
where = 0.577... is the Euler constant. If we account 
for collective effects in the computation of the force auto- 
correlation function, we obtain (see Eq. (96) of Ref. |21|): 



(F(0).F(r)) = (27r) 2 n 



k 2 u{kf 



1 + {2n) 2 f3nu(k) 



e- ikr dk, 



(164) 

where the new terms arise because of the correlations en- 
capsulated in the function h(k). Thus, combining the pre- 
vious results, we get 



(F(0).F(r)) = ^^ +OC 



Jo(kr) 



k 2 



k 2 D 



dk, (165) 



where there is no need to introduce an ad hoc large-scale 
cut-off anymore. We now obtain the result 



(F(0).F(r)> = ^# (M 

Z7T 



(166) 



For kpr — > 0, we recover Eq. ({Ml)) with A = k D \ There- 
fore, the divergence in Eq. I|159(l can be regularized by 
properly accounting for correlations among the particles. 
For comparison, the spatial correlation of the force in 3D 
within the mean-field Debye-Hiickel theory is 



(F(0) • F(r)) = ny 



sin(A:r) 



k 2 



kr 



k D 



dk 



'- ni 2 -e- k °\ 
2 r 



(167) 
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where k]j = (47r/3ne 2 ) 1 / 2 is the Debye shielding length in 
d = 3. Note that the integral remains well-behaved if we 
neglect collective effects (fcc = 0) contrary to the case 
d = 2. On the other hand, the fact that the correlation 
function diverges as 1/r for r — > is due to the fact that 
the distribution W(F) is a particular Levy law, called the 
Holtzmark distribution, whose variance is infinite jHj- 

Let us now consider the temporal correlations of the 
force. If we neglect collective effects, the correlation func- 
tion is given by (see Eq. (93) of Ref. 

(F(0) • F(t)) = (2n) 2 n J k 2 u(k) 2 e~^ vt e~ kH ' '/( 3 « >dk. 

(168) 

Using the same procedure as before, we get 
2 r+°° 



Lt/A 



(F(0) • F(t)) = / ^ZL e -* / W ^ (169) 



which behaves like 



<»<<» ■*<'»-£■■ (5 



(170) 



for wt/yl — > 0. The previous expression indicates that the 
correlation of the force is almost independent of time. This 
implies that the diffusion coefficient, calculated with the 
Kubo formula 

r + 00 r + co 

D= (F(0) • F(t))dt ~ / dt, (171) 
Jo Jo 

diverges linearly with time. This temporal divergence is 
the analogue of the spatial divergence A ~ / + °° d\ ob- 
tained in the kinetic approach of Sec. 14.31 The same ef- 
fect occurs for gravitational systems in d — 3 where the 
logarithmic divergence of the diffusion coefficient can be 
viewed either as a spatial ~ Jdk/k or temporal J dt/t 
divergence (see Sec. 2.9 of HU). Now, if we take into ac- 
count collective effects, the temporal correlation function 
is given by (see Eq. (98) of Ref. 



(F(0).F(i)) = (2^) 2 | 



k 2 u{k) 2 



e(k,k- Vl ) 



re 4k - u */(vi)d Vl dk. 

(172) 



is the time integral of the correlation function (F M {Q)F V (f )) . 
When the dielectric function is taken into account, it is 
seen that the integrals on k in Eq. \ 17311 are now conver- 
gent (see discussion in Sec. 2.8.1 of Ref. |22|)- By these 
means, we can provide a justification of the upper cut-off 
which appears in Eq. I|148|l . The calculations are detailed 
in Appendix IH1 



5 The case d = 1 

For d = 1, the Fokker-Planck equation JJJ reduces to 



dP 
~dt 



= ir(2ir) 2 / dvxdkk^k" 



u{kf 



|e(k,k. Vl )| 



-<5(k.u)/( Vl ), 
(173) 



(27T)'- 



d_ 

dv 



+00 



dv' I dk k 
10 



u{kf 



e(k, kv)\ 2 



S[k(v — v)] 
(174) 



Using S[k(v — v')] = (l/\k\)S(v — v'), the integral over v' 
is straightforward and yields 



— = — A(v) > to,/,- mP^ 



where 



This can be rewittcn in the form 



dP_d_ / dP + p 
dt dv V dv 



with 



and 



D{v)=A{v)Y j m j f j {v) 



rj(v) = —mA(v) 



dfi 
dv 



(175) 



(176) 



(177) 



(178) 



(179) 



By investigating the poles of the integrand (for a Maxwellian 
distribution), it is found that the contribution of each 
mode k should decay with time exponentially rapidly (mod- 
ulated by an oscillatory factor) with an exponent 7^ cor- 
responding to the Landau damping rate. This is the imag- 
inary part of the pulsation u> which cancels out the dielec- 
tric function e(k, to). Therefore, collective effects modify 
the expression of the temporal correlation function and, 
consequently, of the diffusion coefficient. Indeed, accord- 
ing to the Kubo formula, the diffusion coefficient entering 
in the Lenard-Balescu equation 



If the field particles all have the same mass to/, these 
expressions simplify in 

D(v) = A{v)m f f(v), r)(v) = -A(v)mf(v). (180) 

The effective potential Q24H is 

TO 

U(v) = In/(tO, (181) 

TO/ 



and the Fokker-Planck equation can be rewritten 

to d\nf^ 



dP _ d 
dt dv 



D(v) 



with 



D(v) = (2TTYm f f(v) 



(?P_ P _ 

\ dv rrif dv 

50 ku(k) 2 



\e(k,kv)\ 



(182) 



: dk. (183) 
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We note that the stationary solution of the Fokker-Planck 
equation (|182J) is 



P e » = C /"/ («), 



(184) 



and that P eq (v) = f(v) when the mass of the test particle 
is equal to the mass of the field particles m — rrij. More 
explicit expressions of the diffusion coefficient valid for a 
thermal bath are given in |22| . 



6 Conclusion 

In this paper, we have discussed the kinetic theory of sys- 
tems with long-range interactions in a relatively unified 
framework starting from the Landau and Lenard-Balescu 
equations. Using a test particle approach, we have given 
explicit expressions for the diffusion coefficient and friction 
force entering in the Fokker-Planck equation for different 
potentials of interaction in different dimensions of space 
and for different distributions of the bath. We have also 
considered the possibility of a distribution of mass among 
the particles and shown how the results (in particular the 
Einstein relation) are modified in that case. 

For Coulombian and Newtonian potentials in d = 3, we 
have enlightened the connection between results of plasma 
physics and results of stellar dynamics which have been 
derived almost independently and in a relatively different 
form. We have shown that the kinetic equation derived in 
stellar dynamics by Chandrasekhar PI and Rosenbluth et 
al. ^1] from the Fokker-Planck equation can also be ob- 
tained from the Landau 2 equation of plasma physics. We 
have then considered the extension of these kinetic theo- 
ries in d = 2. For a Maxwellian distribution of the bath, 
the diffusion and friction coefficients can be calculated ex- 
plicitly in terms of Bessel functions as shown in [22] • In 
the present paper, we have generalized our approach to an 
arbitrary isotropic distribution of the bath and expressed 
the results in terms of integrals of Bessel functions. In the 
case where the distribution function of the field particles 
is a step function, we have shown that the asymptotic ex- 
pressions of our obtained diffusion and friction coefficients 
reproduce those found by Benedetti et al. [231 ( U P to a fac- 
tor 2(7T— 2)/tt) with a different method. More generally, we 
get a good agreement with their numerical results for all 
velocities. We have shown analytically that, for N — ► +oo, 
the friction diverges logarithmically at the critical veloc- 
ity vo (while for finite N there is no divergence) . We have 
also shown that the linear divergence of the diffusion coef- 
ficient resulting from the Landau approximation could be 
removed by considering collective effects as in the Lenard- 
Balescu treatment of a 3D plasma. Finally, we have shown 
how the results of kinetic theory simplify in d = 1. 

For future perspectives, it could be mentioned that the 
results presented in this paper can be formally extended 
to a generalized class of kinetic equations, associated with 
a generalized thermodynamical framework, introduced by 
Kaniadakis 0Hj and Chavanis |44I39| . These generalized 
equations could be justified in the case of complex systems 
when the transition probabilities from one state to the 



other are different from the form that is usually considered 
due to microscopic constraints ( "hidden constraints" ) that 
affect the dynamics. However, the domains of application 
of these generalized kinetic theories remains to be better 
specified. 



A Alternative derivation of the diffusion 
coefficients and friction force in d = 3 

In Sec. |21 we have derived the expressions of the diffusion 
coefficients and friction force for an isotropic velocity dis- 
tribution of the field particles in d — 3 by using the Rosen- 
bluth potentials. In this Appendix, we show that we can 
obtain the same results by using the method developed in 
Sec. rj] which extends to d = 2. 

Starting from the general expression (|112f> of the dif- 
fusion coefficient and introducing a spherical system of 
coordinates with the z-axis in the direction of v, we find 
after straightforward algebra, that the diffusion tensor can 
be written as in Eq. I|18|l with 



D\\ = (2 7 r) 7 Vm 1 1 



-\-oo 



dx 



2 1 



3 

2 \ sinx cos a; 

■ 4 — =- 



k 3 u(k) 2 dk 



(;) • < 185 > 



1 f +oc 

D± = (2tt) 7 Vm ? - / k 3 u(k) 2 dk 
, Wo 



sin x cos x 
4 — = 4 — r- 



(186) 



The Fourier transform of an isotropic velocity distribution 
of the field particles in d = 3 is 



0(k) = 



1 



27T 2 fc 



f{vi) sm(kvi)v\dvi 



(187) 



Substituting this expression in Eqs. i|185|) and l|18fci|) and 
introducing the functions 



dx 



2 \ siM cos a; 
2 1- -=■ +4- 



— sin(Air), 
x 

(188) 



r+oo 

I±(X)= I d.r 



sin x cos x 
4 — = 4- 



■ sin(Air), 



(189) 



we obtain 



D|, = 8A Vm, f fj(vi)I\\ (—) vxdvx. (190) 

3 V 

D ± = SAj^m, f + fjiv^Ix v x dv x . (191) 

T « ^ 
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The integrals ((188(1 and 1(189(1 can be calculated explicitly. 
Using the results: I\\ (A) = vr/3, Jj_(A) = 2tt/3 for A > 1 
and 7||(A) = (tt/3)A 3 , 7j_(A) = ttA - (tt/3)A 3 for A < 1, we 
recover the expressions ((86(1 and 1)87(1 . The friction r\ can 
be obtained from Eq. 1(20(1 and we recover Eq. 1(89(1 . 



B Asymptotic behaviors of the functions 
dHZb and (HlHl) 

B.l Asymptotic behaviors for A — > 

In this Appendix, we determine the asymptotic behaviors 
of the functions 1(127(1 and ((128(1 for A — > 0. Using identities 
ifTTHl and jnU), Eq. lfHT)) can be rewritten 



R n (\) = - f decode [ °° e lxcos9 ^^dx. (192) 
^Jo & 



Setting i = cos 6*, we obtain 

fi„(A) = i r 



+1 t 2 dt r +aa 



e'^fc (193) 



In this expression, t and x are real and the domains of 
integration tq : — 1 < t < 1 and Co • < ^ < +oo are 
taken along the real axis. Under these circumstances, we 
cannot simply expand the last function in Taylor series for 
A — > because the resulting integrals would not be con- 
vergent. However, regarding x and t as complex variables, 
it is possible to choose paths of integration along which 
this expansion will converge. We shall first carry out the 
integration on x, for a fixed t. It will therefore be possible 
to choose the (complex) integration paths for x dependent 
on t. The integration paths are modified as follows: tq is 
replaced by r, the semi-circle with radius unity lying in 
the domain T m (t) > 0. Therefore, arg(i) varies from it to 
when t moves from —1 to +1. On the other hand, Co 
is replaced by C^t j the line starting from the origin and 
forming an angle 



with the real axis. With these new contours, 



(194) 



1 f t 2 dt 
En (A) = -Re / - 



e i Xt J±(M dx (195) 



When t moves from —1 to +1, tp t varies from — ^ to +~. 
The argument of ixt is equal to 7r. Therefore, the real 
part of ixt is always negative and the function e lxt decays 
exponentially to zero as \x\ — > oo. Therefore, with the new 
paths of integration r and Ci/> t , it is possible to expand the 
integrand of Eq. 1(195(1 in power series of A, for A — > 0, and 
integrate term by term: 



a, (A) 



-Re 



t 2 dt 



3^,2 



X s x 

IF 



dx. 
(196) 



Setting ixt = —y (y real > 0), we get 



#H (A) 



1 



-Re 



-1 + 2 



t 2 dt 



x / e 
/o 



7T 7-i VT^T 2 

2 16t 2 it' 



(197) 



where we recall that t is a complex variable and the inte- 
gration has to be performed over the semi-circle of radius 
unity lying on the domain I m {t) > 0. Writing t — e l9 , we 
find that 



and 



Therefore, 



f +1 tdt 
Re i / ■ = 0, 



Re 2 / ; = 7T. 



Using a similar procedure, we find that 
fi ± (A) ~ ^. 



(198) 



(199) 



(200) 



(201) 



Note that this second result can also be obtained directly 
from Eq. 1(128(1 by using J\(\x)/x ~ A/2 and 



Ji(ar) 



rfa; = 1. 



(202) 



B.2 Asymptotic behaviors for A — > +oo 



We now determine the asymptotic behaviors of the func- 
tions ((127(1 and l(128|) for A — > +oo. Setting z = Ax, we 
can rewrite Eq. 1(127(1 as 



Using identity ((114(1 . we obtain 



(203) 



k Jo 



d9 sin 2 6 



+ 00 







(204) 



Setting < = cos#, using the contours introduced in the 
previous section and expanding the last term in Taylor 
series for 1/ A — * 0, we find that 



J?H (A) = -Re / dty/l - t 2 [ 



1 _ 3_z^ 

2 _ 16 X 2 



dz. 
(205) 
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Setting izt = —y (y real > 0) we get 



We also note the particular values 



J?,, (A) 



1 

—Re 



dWi-t 2 



I e 

'o 



1 3 y 2 

2 + 16 \H 2 



% (206) 
it 



where the integration on t has to be carried out on the 
upper semi-circle in the complex plane. Using the identi- 
ties 



Re i 



dt 



I —UL 

VI - t 2 — = 7T 



Re ; / \J\-t 2 13 



dt 7T 

_ ~ ~2' 



wc find that 



Using a similar procedure, we obtain 
fi±(A) = 1 ' 



16A 2 



(207) 



(208) 



(209) 



(210) 



B.3 Relation to hypergeometric functions and 
behavior for A ~ 1 

In fact, the functions H127|l and 1)128(1 can be expressed 
in terms of Hypergeometric functions. We can then easily 
obtain their asymptotic behaviors from standard formulae 
|35) . We have 



i?H (A) 



*x(A) = ^(-±,i,2,A3 



For A -> 0, 



fill (A) 



A 3 A 5 



15A 7 



32 1024 



R±(X) = — — — — — 

K ' 2 16 128 2048 



and for A 



fill (A) 



16A 2 128A 4 



1 1 



1 



2 16A 2 128A 4 



(211) 
(212) 

(213) 
(214) 

(215) 
(216) 



3tt 

From Eq. I|211|) . we get 



fi||(l) = — , R x (l) = —. 



3tt 



(217) 



fill (A) 



A 2 " 



(218) 



This function diverges for A — * 1 like 



1 



fill (A) = (2 7 + ln2 + ln|l- A| + 2^(3/2)), 



37T 7T 



(219) 



where 7 = 0.577216... is the Euler constant and ip(3/2) = 
0.03649... is the Digamma function |45|. 



C Regularization of the linear divergence 

In this Appendix, we show how the linear divergence of 
the diffusion coefficient for a Coulombian plasma in d = 2 
can be regularized by taking into account collective effects. 
When collective effects arc taken into account using Eq. 
instead of Eq. JSJ, the diffusion coefficient is given by 



D"" = ir(2ir) d m / dv 1 dkk^k 1 ' 



u{k) 2 



|e(k,k-v)| 2 



5(k-u)/(vi). 

(220) 

We concentrate here on a thermal bath with Maxwellian 
distribution. Using the same method as in Sec. 12.51 but 
keeping the dielectric function, we obtain 



= 7r(2n) d mp ( ^ 
2n 



1/2 



x / dk 



k |e(k,k-v)| 2 
For a thermal bath one has (see, e.g., \22\): 

|e(k, k ■ v)| 2 = (1 - V (k)B(k ■ x)) 2 + C(k • x) 2 , 



(221) 



(222) 



where we have defined k = k/fc, x = (/3m/2) 1 / 2 v, <q{k) 
-{2ir) d u{k)f3mp, B{z) = 1 - 2ze~ z2 /* e e 'dt and C(z) 
y/Tt\z\e~ z . The diffusion tensor can be rewritten 



= 
where 



{2ii) d j3 2 pm V 2tt 



(223) 



1/2 



I(z) 



+ OO 



k d j](k) 2 dk 



(l-r/(fc)B(z)) 2 +C(z) : 



(224) 
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For a Coulombian potential, one has rj(k) 
we obtain 



-k 2 D jk 2 and 



I(z) 



k d dk 



(B(z) 



k 2 /k 2 D ) 2 



C(z) 



(225) 



This can be rewritten I(z) — k'^ 1 C(z) L 
where $ d (z) = / + °° t d dt/((z + t 2 ) 2 ' 



collective effects we have instead 



r$ d ( B (z)/C(z)) 
1). If wc neglect 



■idau 



dk 



(226) 



which presents divergences for k — ► if d < 3. The reg- 
ularization of the divergence in d = 3 has been treated 
by Balescu [8] and the regularization in d = 1 has been 
treated in [22]. Let us focus here on the case d = 2. We 
note that, contrary to Eq. I|226|l . the integral (|225|l is well- 
behaved for k — > and k — > +00. Therefore, there is no 
linear divergence of the diffusion coefficient and friction 
force when collective effects are taken into account. 

To obtain the expression of the diffusion tensor, one 
has to substitute Eq. (|225|l in Eq. (|223|l . using z = xcos8, 
and carry out the integrations. We shall simplify the cal- 
culations a little bit by introducing an approximate an- 
alytical expression of I(z) (we have checked numerically 
that the exact treatment yields close results). Let us first 
consider asymptotic behaviors of the previously defined 
functions. We have <P 2 (z) = (tt/(2V2))(1-z/2 + z 2 /8 + ...) 
for z — > 0, #2(2) ~ 7r/(4yi) for z — > +00 and # 2 (z) ~ 
(ir/2)\/—z for z — » —00. On the other hand, B(z)/C(z) ~ 

l/(v^r|*|) for z ^ and B(z)/C{z) l/(20F|z| 3 )e 22 

for |z| — > +00. Thus, we find that I(z)/k 3 D — > 7r/4 for 
z -> and I(z)/k 3 D -> 0r/(2V2z 2 )e z2 for 2 -> +00. Let 
us consider a simple interpolation formula of the form 

J(z) 



1,3 



7T 

4 



1/2 



(227) 



With this expression, it turns out that the trace of the 
diffusion tensor 



=D I(xcos6)e- x2 
Jo 



9 d6, 



(228) 



can be calculated analytically (Do is the value of the con- 
stant in front of the integral in Eq. H223(l L When we use 
Eq. H227[) . we obtain 



IT 3 / 2 2 



(229) 



Alternatively, when we use Eq. 



we get 



Landau 



= D k%2irAe-^ I [ ~— ) , (230) 



where A — J Q °° dk/k 2 . We note that the trace of the 
diffusion coefficient decays like x" 1 in each case. We can 






, . 1 


\ D Landau (x)/D k D 4 A 




, \ 




V D(x)/D k D 3 






k D A(x) 







10 



Fig. 8. Dependence of the trace of the diffusion tensor with 
the velocity, using (i) the Landau approximation or (ii) the 
Lenard-Balescu treatment of collective effects. 



write Eq. 



in the form of Eq. it53U|) : 



DW = D k%2nA(x)e--I — 



(231) 



where A(x) is now a function of x which is perfectly well 
defined (without divergence). Using Eq. I|229l) and com- 
paring with Eq. (|231|l we find for x — > that 



and for x 



-00 that 



-A(+oo) 



/27r 1.412. 



4fc D 



ko 



(232) 



(233) 



This justifies, without introducing an ad hoc large-scale 
cut-off, that A is of order the Debye length k^ 1 as ex- 
pected. We also note that the function 



k D A(x) = - 



h(x 2 /2) 
h{x 2 /2) 



(234) 



does not vary crucially with x and remains of order unity 
(see Fig. IHJ. Therefore, using the Landau approximation 
and introducing a large scale cut-off at the Debye length 
kp 1 seems to be a reasonably good approximation. How- 
ever, when one evaluates the component D»(x) of the dif- 
fusion coefficient, one finds that the treatment using Eq. 
(I225|) leads to a decay like x~~ 2 for x — > +00 while the 
Landau approximation i|22ti|) leads to a decay like a; -3 
(see Sec. 14.20 . Therefore, in that case, there are qualita- 
tive discrepencies between the two approaches. 
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